Correlated-Photon Imaging with Cancellation of Object-Induced Aberration 
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We show that a recently discussed apparatus for aberration-cancelled interferometry may be 
modified to perform correlated-photon imaging in the so-called "ghost" imaging configuration. For 
objects in the vicinity of a particular plane, the images are free of object-induced phase distortions. 
This apparatus has the distinctive feature that it may be used to superimpose images of two objects 
in a manner that could lead to useful effects and applications. We show that the apparatus works 
using either quantum-entangled or classically correlated light sources. 

PACS numbers: 42.30.Va,42.15.Fr,42.30.Kq 



I. INTRODUCTION 



Correlated-photon imaging, sometimes known as 
"ghost" imaging, was first discovered using entangled 
photon pairs [l[ from spontaneous parametric downcon- 
version (SPDC). It has since been found that most as- 
pects of ghost imaging can be simulated using spatially- 
correlated classical light 0, including thermal and 
speckle sources A separate line of research has 

shown that entangled photon pairs from downconversion 
may also be used to cancel some of the effects of frequency 
dispersion [Tol - [T2j or spatial dispersion (aberration) 
IToT ]. In [lU, it was pointed out that it is possible to con- 
struct an interferometer such that if an object is placed 
in a particular plane then the effects of all phase shifts 
induced by that object, including all object-induced aber- 
rations, will cancel in the resulting coincidence rate. The 
goal here is to move away from interferometry and to 
produce an analogous effect in an imaging system. We 
show that this may be achieved by a simple variation of 
the traditional ghost imaging apparatus of figure 1. It 
is thus possible to produce images of the object's ampli- 
tude transmittance profile, undistorted by phase effects 
as long as the object is entirely contained within a small 
region near the special plane mentioned above. (For sim- 
plicity we will only discuss transmission here; the case of 
reflection at the object is similar.) We then show that, al- 
though an entangled source was required for the temporal 
correlation experiments discussed in [l3rll5j . a classical 
source with transverse spatial correlation will suffice for 
imaging. 

In addition, if two object are placed in the resulting 
optical system, one in each arm, the image produced will 
simply be the point-by-point product of the images that 
would be generated by each of the two separately. This 
is a new feature that does not appear if two objects are 
placed in the arms of other types of ghost imaging sys- 
tems. We will comment on several possible applications 
of this effect below. 

We begin by briefly reviewing ghost imaging in section 
2, followed by a review of aberration-cancelled interfer- 
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FIG. 1: (Color online) Schematic depiction of correlated- 
photon imaging setup. The photons in the two arms have 
anticorrelated transverse momenta ±q. 



ometry in section 3. We then show how a small change 
converts the aberration-cancelled interferometer into a 
new type of ghost imaging system. We analyze this imag- 
ing system first with an entangled light source in section 
4, then with a classical source in section 5. Finally, in 
section 6 we discuss an important technical point about 
the need for lenses in front of the detectors, followed by 
conclusions in section 7. 



II. CORRELATED-PHOTON IMAGING 

Correlated-photon imaging or "ghost" imaging |l| is 
done with an apparatus like the one depicted schemati- 
cally in figure 1. In the original version, the correlated 
photon source is a nonlinear crystal pumped by a 
laser, leading to SPDC. Entangled photon pairs with an- 
ticorrelated momentum components q and — q transverse 
to the propagation direction travel along the two arms of 
the apparatus. The object to be viewed is placed in arm 
2 (the upper branch), followed by a bucket detector, D^. 
D2 can not record any information on the position or 
momentum of the photon that reached the object; it can 
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only tell us whether or not the photon reached the de- 
tector unimpeded by an object. The other arm has no 
object, and all the photons reach a CCD camera or ar- 
ray of pointlike detectors without hindrance. A lens may 
be inserted in this branch for image formation. A co- 
incidence circuit is used to record a count every time a 
photon detection occurs simultaneously (within a short 
time window) at each detector. By plotting the coin- 
cidence rate as a function of position Xi in detector 1, 
we build up an image of the object. This is true even 
though photons that actually encountered the object in 
branch 2 left no record of the object's position, and the 
photons in branch 1 that do carry position information 
never encounter the object. 

The crucial ingredient is the spatial correlation of the 
photon pair. It was found @, Q that entanglement was 
unnecessary: a classical source with anticorrelated trans- 
verse momenta could mimic the effect. The correlated 
light source in this case consists of a beam steering mod- 
ulator (a rotating mirror, for example) directing a classi- 
cal light beam through a range of q vectors, illuminating 
different spots on the object. The beamsplitter turns the 
single beam of transverse momentum q into a pair of 
beams with momenta q and — q. The results were simi- 
lar to those with the entangled source, but with half the 
visibility. It was later shown that thermal_and speckle 
sources may also lead to ghost imaging 



III. SUMMARY OF ABERRATION 
CANCELLATION IN QUANTUM 
INTERFEROMETRY 



Consider the setup shown in figure 2 [13Hl5l | . Each 
branch contains a 4f imaging system with lenses of focal 
length / and a thin object that provides spatial modula- 
tion Gj(y) of the beam, where j = 1,2 labels the branch, 
and y is the position in the plane transverse to the axis. 
The goal is to cancel object-induced optical aberrations 
(position-dependent phase shifts produced by the Gj). 
The case of a single object in one branch only may be in- 
cluded by simply setting G = 1 in the other branch. The 
plane of the samples (labelled II in fig. 2) is the Fourier 
plane of the 4f system. Time delay r is inserted in one 
branch. In the detection stage, two large bucket detec- 
tors Di and D2, connected in coincidence, record the 
arrival of photons, but not their positions. Apertures de- 
scribed by pupil functions pi(xi) andp2(x2) are followed 
by crossed polarizers at 45° to each beam's polarization, 
before arriving at the detectors. 

A continuous wave laser pumps a nonlinear crys- 
tal, leading to collinear type II parametric downconver- 
sion. The frequencies of the two photons are firj ± v 1 
with transverse momenta ±q. For simplicity, assume the 
frequency bandwidth is narrow compared to firj- The 
two photons have total wavenumbers a ° ±v k, The 
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FIG. 2: (Color online) Schematic of interferometer with even- 
order aberration-cancellation. Large bucket detectors Di and 
D2 are integrated over and connected by a coincidence circuit. 



downconversion spectrum is 



■ J,A(q,y) 



where L is the thickness of the crystal, and 



A(q, v) = -i/D + Me 2 -q- 



2|q| : 



(1) 



(2) 



vpump 



D is the difference between the inverse group velocities of 
the ordinary and extraordinary waves in the crystal, and 
M is the spatial walk-off in the direction e"2 perpendicular 
to the interferometer plane. The last term in A is due 
to diffraction in the crystal. Ignoring the vacuum term 
and terms of higher photon number, the wavefunction 
entering the apparatus is approximately given by 

I*} = J d 2 q dv $(q,i/) 

x at( q} n + i/)at(-q l n -i')|o) > (3) 

where a s and hi are annihilation operators for the signal 
and idler photons. For collinear pairs, horizontally polar- 
ized photons are directed into the upper branch and ver- 
tically polarized photons into the lower branch by means 
of a polarizing beamsplitter. Alternatively, noncollinear 
pairs could be used with polarizers selecting horizontal 
(H) polarization in the upper branch and vertical (V) in 
the lower one. In either case, we will refer to the H pho- 
ton in the upper branch (branch 2) as the signal and the 
V photon in branch 1 as the idler. 

The coincidence rate is of the generic form [l6[ 



R(r)=Ro l-AM-ig, 
where A(x) is the triangular function 



W(t) 



l-\x\, 
0, 



\x\<l 
x\ > 1 



(4) 



(5) 
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For large apertures, pi(xi) = ^2(^2) ~ 1; so, as shown 
in [14| . the background and r-modulation terms are 



Ro 
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(6) 
(7) 
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where k is the longitudinal wavenumber. 

We may write Gj(x) = tj(x)e i ^^ x ', with tj real. 
Aberration effects arise from spatially-dependent phase 
factors <pj(x), which lead to distortion of the outgoing 
wavefronts. The phase functions may be decomposed 
into a sum of pieces that are either even under reflec- 
tion, ^ en) (-x) = ^ en) (x) or odd, cj>f dd \-x) = 

—<jy° {x). Astigmatism and spherical aberration, for 
example, are included in the even-order part, whereas 
coma is odd. 

In equation Q, the factors G\ Gi (-^) be- 



(8) 



the form of the difference in the exponent shows that 
even-order aberrations arising from object 1 cancel from 
the modulation term. The even-order aberrations from 
object 2 cancel in a similar manner. This is the even- 
order cancellation effect demonstrated in 13 1 and [i~4j . 
As pointed out in fl5l ). even- and odd-orders cancel si- 
multaneously only in the special case G\ = G%. These 
cancellations are exact only for aberrations induced by 
thin objects in the particular plane EL 

In the term i?o, both even-order and odd-order aber- 
rations cancel even for G\ 7^ G2. For time correlation 
experiments, this is an unimportant background term; 
however this term is the foundation of the imaging ap- 
paratus described below, since the beamsplitter will be 
absent, meaning that there will be no modulation term 
W(t). The physical mechanism of the various possible 
cancellations are discussed in more detail in 15 1. 



IV. ABERRATION-CANCELLED GHOST 
IMAGING WITH ENTANGLED SOURCE 

Now we wish to look at the ghost imaging analog of 
the aberration-cancelling interferometer of the previous 
section. This leads us to the hybrid device of fig. 3. 
This new apparatus differs from that of fig. 2 in several 
respects. First, we have removed the time delay, polar- 
ization filters, and beam splitter; these were needed to 
produce the interference effects desired in [l3Hl5l ]. but 



FIG. 3: (Color online) Schematic of correlated-photon imag- 
ing setup with aberration- cancellation. All orders of aberra- 
tion cancel. 



are not necessary for imaging purposes. Also, in order 
to obtain spatial resolution, one bucket detector (D%) 
is replaced by a moveable pointlike detector or a CCD 
camera. The removal of the beam splitter and the intro- 
duction of spatial resolution are the key changes. After 
allowing for an arbitrary source of correlated (quantum 
or classical) light, we arrive at an apparatus in fig. 3 
that looks very much like the ghost imaging setup of fig. 
1, but with a 4f imaging system in each branch. In this 
section, we assume that the light source is parametric 
downconversion. 

The coincidence rate at location Xi of D\ is 



-ff(xi) = J d 2 X2dt 1 dt 2 \A(x 1 ,x 2 ,t 1 ,t2) 



(9) 



where the transition amplitude is 

A(xi,x 2 ,t 1 ,t 2 ) = {0|^ +) (x 1 ,t 1 )^ +) (x 2 ,t 2 )|*). (10) 

Taking the two detection apertures described by p\ and 
P2 to be large, we compute the coincidence rate to be 



i?( Xl ) = {[B(x 1 ) + B(-x 1 )] + [G(x 1 ) + C*(x 1 )]} 
... / / \ ,-. / ./ \ 

x 



where 



fcxi \ / fcxi 



(11) 

(12) 
(13) 



Using eq. ([T]), these integrals may be evaluated; they 
turn out to be xi-indepcndent constants. Sweeping all 
overall constants into a single constant TZq, we find: 



(14) 
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Only the modulus of each Gj enters into i?(x), so we 
see that the aberrations introduced by the object phases 
cancel to all orders. This will be exact only in the Fourier 
plane, but as was true for the interferometer case, we 
would expect it to continue to remain approximately true 
as we move out of the plane up to a maximum distance 
on the order of where / and a are the focal length 
and radius of the lens and r s is the maximum radius of 
the object being viewed. (See ref.([l5j]) for a derivation 
of this estimate.) 

If Gi = 1, then we have an ordinary (non-ghost) image 
of |Gi|. On the other hand, if G\ = 1 then we have an 
inverted ghost image of |C?2 1 ■ In either case, the image is 
magnified by a factor of m = fy-. Note that, in contrast 
to the interferometry case, even and odd order phases 
both cancel, even in the general case G\ ^ Gi. The can- 
cellation of the phases only occurs in the Fourier plane; 
aberrative distortions begin growing when the objects are 
moved out of this plane. 



V. APPLICATIONS TO IMAGE ANALYSIS 

Note that if both G\ and Gi are nontrivial objects, 
what we actually see is their pointwise product. This is 
a distinctive feature of this apparatus. It can be verified 
by straightforward calculation that the simple product 
structure of eq. (1141) does not occur in other obvious vari- 
ations of two-object ghost imaging systems; for example, 
it does not occur if the 4f imaging system in either branch 
(or both) is replaced by a single lens imaging system, or 
in a system without lenses. We may make use of this 
product structure in a number of ways. For example, if 
the object of interest is G2 but the second branch intro- 
duces some known distortion to its image (e.g. there may 
be aberration in the optical system or variations in the 
refractive index of the propagation medium), then this 
can be cancelled by using an object G\ that introduces 
an opposite distortion (via a deformable mirror, for ex- 
ample). The two distortions then cancel, leaving no net 
effect on the image. Alternatively, if G2 has a dim, low- 
transmissivity area that we wish to view, but it is being 
obscured by a bright, high-transmissivity area nearby, we 
may use a mask for G\ which allows a view only of the 
twins of photons coming from the dim region of inter- 
est, blocking photons that are partnered with light from 
other areas. 

At first glance, the product structure might seem to 
open up a further interesting possibility. Suppose G\ 
is the object we wish to view. If there is no object in 
branch 2 (G2 is simply a constant), then the resolution 
with which we may view G\ is the same as if the sec- 
ond branch was not there. It would be limited by the 
sizes of the Airy disks produced by the lenses. However, 
if G2 is taken to be a small pinhole, the area we would 
be able to see of G\ at any given time would be lim- 
ited by the size of the pinhole. Thus, it would seem that 
by taking the pinhole small enough, we would be able to 



limit our view of Gi to an area smaller than the standard 
Abbe limit, thus achieving subresolution imaging. Unfor- 
tunately, when the finite sizes of the lenses are properly 
taken into account (eq. (fT4"]l was derived in the limit of 
large lenses), the combined action of diffraction in the 
two branches conspires to give the single-branch resolu- 
tion as its best-case limit, occurring when the pinhole 
radius is negligible. As the pinhole radius at G2 grows 
to finite size, the resolution becomes worse than in the 
single-branch case. 

One additional observation on applications of the prod- 
uct structure arises if we replace the position-resolving 
detector in branch 1 by a bucket detector, thus introduc- 
ing an integration over xi . We have now lost all imaging 
ability, but note what happens if we displace one of the 
objects (object 1, say) by some distance in the transverse 
plane. If the 2-dimensional displacement vector is r, then 
equation (TJU) is replaced by 



Gi 



J_ 

fn 



(xi + r) G 2 - 
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fi 



-xi 



dx\. 



Thus, despite 
tor has spatial 
cally computes 

ff( r ) = (h 

the magnification. 



(15) 
detec- 
opti- 



the fact that neither 

resolution, the system 
the spatial intensity correlator 
(x + r)) I 2 (— m _1 x))}, where m is 
(The correlation here is actually 
between the object Gi(x) and the inverted object 
G2(— x), but an additional lens can be added to remove 
remove the inversion and cancel the minus sign in G2.) 
The full correlation function can be found by moving one 
object repeatedly to scan over the full range of relevant 
r vectors. Taking one of the two objects to be unknown 
and the other to be some known template, this could 
provide a means of identifying the unknown object by 
quantifying its degree of similarity to the template. This 
could be useful, for example, in comparing silicon chips 
on an assembly line to a standard chip, and identifying 
those chips with flaws. Note in particular, that the un- 
known object may be in a remote, inaccessible location; 
for example, the object might be a cell inside the body 
being viewed through an endoscope and compared to a 
cell in the lab. As in the case of the temporal correlator 
studied with the interferometer of refs. EMU, the 
effect of object-induced aberrations (differences between 
phase shifts induced by the two samples) cancels out of 
the spatial correlator. 



VI. IMAGING WITH A CLASSICAL SOURCE 

We now replace the downconversion source of the pre- 
vious section by a classical source of anticorrelated pho- 
tons, as in [2]. Light entering a beam splitter with trans- 
verse momentum q leads to outgoing beams with anticor- 
related momenta q and — q. If the beam steering modula- 
tor produces momentum spectrum /(q), the input state 
for pairs of photons having the same q before the beam- 
splitter will be ~ J cP<7F(q)a£(q)a],(q)|0}, where a£ is the 
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creation operator for pump photons and F(q) = / 2 (q). 
We assume for simplicity that -F(q) is an even function, 
F(q) = F(— q). Denoting creation operators in the two 

outgoing branches by a\ and a\, the incoming photon 
pair will produce a state after the beamsplitter given by: 

|*> = \J <*V(q) [aI(q) + o|(-q)] (16) 

X [ot(q) + aJ(-q)] |0) 

= J dV(q) [a|(q)a+(-q) + . . . ] ]0), (17) 

where the terms dropped in the last line are those which 
do not contribute to coincidence detection. The detection 
amplitude of eq. (|10|) is then proportional to 



d 2 qF(q)e 



iq-(x!-x 2 ) 



ffi(q,Xi)ff 2 (-q,X2), (18) 



where Hj is the transfer function for branch j. Integrating 
over D2, we then have the coincidence rate: 



-R(xi 



I -7-x, ) r;, 



xi G 2 



Id 



Xl 



(19) 



This is similar to the result for the entangled-source appa- 
ratus, except modulated by the factor F (j^ x lj which is 
determined by the details of the beam steering modulator. 
Similarly, for thermal or speckle sources, this factor will arise 
from the transverse momentum spectrum of the source. 



VII. ROLE OF THE DETECTION LENS 

Consider now the lenses immediately before the detectors 
in fig. 3. With no such detection lens present, the transfer 
function for branch j would be 



k 



(20) 



from which we see that the information from each q value is 
spread over all x values. But with the lens, eq. (|20[) becomes 



Hjfaxj) = e m W Ve" 



• <■ (%W^ <->j.) 



k J \fo 



so that each q value is localized at a single point in the de- 
tector plane via the delta function. Since each q value is also 
matched to an object point, the localization in the second 
case defines a mapping between points in the object plane 
and points in the detection plane, allowing reconstruction of 
an image by the pointlike detector D\. This can be verified 
by computing the coincidence rate with or without the final 
lenses, i.e. using either eq. (|20[) or eq. (|21[) . Doing so, we find 
that without the branch 1 lens the coincidence rate becomes 
independent of xi, making imaging impossible. In contrast, 
removing the branch 2 lens has no effect. This makes intuitive 
sense: we integrate over x 2 , so it does not matter if the mo- 
mentum information in this branch is localized or spread over 
the entire detector. Thus we arrive at an important technical 
point: the lens before the bucket detector may be removed 
without harm, but the branch 1 lens is essential for imaging. 

The need for a lens before D\ may be viewed as follows. 
The 4f system in either branch transfers modulation Gj from 
the transverse coordinate space (x) to the Fourier space (q), 
which is where the aberration cancellation actually takes place 
(see [H|]). The lens in front of D\ is then needed to transfer 
the modulation back to coordinate space for imaging. 



VIII. CONCLUSIONS 

In conclusion, we have proposed a new type of two-object 
ghost imaging apparatus that cancels phase effects from thin 
objects in the vicinity of a particular plane and which allows 
comparisons between pairs of objects. The method involves 
a relatively simple apparatus and can be done with either 
entangled photon pairs or with a classically-correlated light 
source. This apparatus may have potential for new applica- 
tions in biomedical research, industry, and other fields. 
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